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The linear perturbation group transformation (LPRG) is used to study the thermodynamics of the axial next- 
nearest-neighbor Ising model with four spin interactions (extended ANNNI) in a field. The LPRG for weakly in- 
teracting Ising chains is presented. The method is used to study finite field para-ferrimagnetic phase transitions 
observed in layered uranium compounds, UAsi-j^Se^, UPd2Si2 or UNi2Si2. The above-mentioned systems are 
made of ferromagnetic layers and the spins from the nearest-neighbor and next-nearest-neighbor layers are 
coupled by the antiferromagnetic interactions /i < and J2 < 0, respectively. Each of these systems exhibits a 
triple point in which two ordered phases (ferrimagnetic and incommensurate) meet the paramagnetic one, and 
all undergo the high field phase transition from para- to ferrimagnetic (-i- -i- -) phase. However, if in \iksi-xS&x 
the para-ferri phase transition is of the first order as expected from the symmetry reason, in UT2Si2 (T = Pd, Ni) 
this transition seems to be a continuous one, at least in the vicinity of the multicritical point. Within the MFA, the 
critical character of the finite field para-ferrimagnetic transition at least at one isolated point can be described 
by the ANNNI model supplemented by an additional, e.g., four-spin interaction. However, in LPRG approxima- 
tion for the ratio )c = /2//1 around 0.5 there is a critical value of the field for which an isolated critical point 
also exists in the original ANNNI model. The positive four-spin interaction shifts the critical point towards higher 
fields and changes the shape of the specific heat curve. In the latter case for the fields small enough, the specific 
heat exhibits two-peak structure in the paramagnetic phase. 
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1. Introduction 

The Linear Perturbation Renormalization Group (LPRG) (l| method uses a simple one-dimensional 
decimation to study universal (critical) and non-universal such as a location of the critical temperature 
and temperature or field dependence of the properties of thermodynamic quantities of several classical 
and quantum higher-dimensional models. For the first time, this kind of method was proposed by Suzuki 
and Takano (ST) |2]. In the ST approach the one-dimensional decimation is combined with the Migdal- 
Kadanoff (MK) bond moving approximation. The disadvantages of the latter method especially for the 
quantum systems were discussed by Barma et al. [St] and Castellani et al. |4J. Here, we wish only to remind 
that the MK approach gives rather poor quantitative results even for the two-dimensional Ising model 
and there is no possibility to construct any systematic approximation procedure within this method. 
For example, the MK procedure gives for the Ising model on the square lattice the values of the inverse 
critical temperature kc ~ 0.61 whereas the exact value is kc ~ 0.44 and for the s = | XY model kc ~ 
1.2 01 much larger than the value kc ~ 0.64 estimated from the high-temperature series expansion 0], 
kc ~ 0.71(67) found from Monte Carlo simulations by fitting to the exponential law (^[^and power law 
@l, respectively or rotationally invariant non-linear (block) transformation fcc ~ 0.62 The LPRG 

method has been proposed for the so-caUed quasi-one-dimensional magnets made of spin chains with 
the intrachain coupling k and much weaker interchain coupling ki < k. However, even for the standard 
Ising model fci = k, the LPRG approximation gives the results which are in very good agreement with the 
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exact ones kc ~ 0.45 . For fc > fci > O.lSfc, the deviation from the critical temperature exact values is less 
than 2% and for 0.5A; > fci > 0.15fc even less than 1% [11.1 . 

In this paper, the LPRG is used to study the thermodynamics and the existence of a critical point in 
two-dimensional axial next-nearest-neighbour Ising model with four spin interactions (extended ANNNI) 
in a field. 



2. LPRG 

We first describe in detail the LPRG approach in the simplest possible case, i.e., the Ising chains with 
intrachain interaction / coupled by the weak interchain interactions /i defined by the Hamiltonian 

^ = fc ^ SijSij+i + fci X SijSi+ij , (2.1) 
('•;■> ('7> 

where the label i refers to rows and j refers to columns, the factor - l/k^T has already been absorbed 
in the Hamiltonian (fc = J/k^T, fci = Jilk^T), and fci < fc. The LPRG approach starts with an exact deci- 
mation for one-dimensional system. Dividing the spins of the chain into two groups, in the simplest case, 
S2i+\ (odd spins — survived) and S2i+z (even spin — decimated) one can write the decimation transfor- 
mation in the form 

g^(S2i + l) _ -pj.^ g^(S2i + i,S2i+2) (2 2) 

In each step of the transformation d2.2> . every other spin is decimated. The same transformation can be 
written by using a linear weight operator P((t,-, s;) 

e-^''''^ =TrsP(cr,s)e'^''\ (2.3) 

where the weight operator P{a, s) which couples the original (s) to the new spins [a) is chosen in a linear 
form as follows: 

1 ^ 

P{(J,s)^ — Y\{\ + Oi+iS2ui). (2.4) 

The weight operator P{o,s) projects the original spin s space onto the space of the effective spins 
o. For the one dimensional Ising model with interchain interaction fci - 0, the transformation \2.2\ or 
equivalently <2.3> can be carried out exactly. So, we can separate the Hamiltonian 1 12. 1> in a manageably 
exactly unperturbated part J6'q containing intrachain interaction fc [first term of the Hamiltonian I I2.1H , 
and a remainder ^/ containing interchain interaction fci [second term in equation <2.1) 1. With the nota- 
tion 

zo =TrsP(cr,s)e'^''''' (2.5) 

and 

<^)o = —TisAPia, sje-^"'"' (2.6) 



the transformation <2.3) can be written as 

^'(cr) =lnzo + ln<e'^''">, (2.7) 
with the following cumulant expansion for ln<e'^''*'> Jlill 

ln(e^/W> = {jei)Q + ((^/)o - <^/)g) . (2.8) 

The idea of LPRG in two dimensions is presented in figure [l] The chains are divided into two black 
and white groups. In each renormalization step, every other spin from a black chain is decimated and in 
a white chain aU spins are removed. As a result, one gets the system of effective spins o. The single chain 
"partial" partition function zq can be easily found as 11311 

zo - (cosh2fc-i- 1) -I- (cosh2fc- l)(7ijCT/j+i , (2.9) 



13704-2 



LPRG for Ising-like systems 



0-1,1 0■^,2 



Sl,1 Sl,2 Sl,3 

• • • 

o o o 

• • • ■ ■ 

o o o 

Figure 1. The LPRG procedure for weakly interacting cliains in two dimensions. 

and 

1 1 
Inzo = -(ln2 + lncosh2fc) + — cr ij a ij +i\ncosh2k. (2.10) 



To evaluate the cumulants <2.8K one has to know the averages (sij) and (stj . . . sij+n) from the black 
(decimated) and white (removed) rows. For the spins from the black rows (figure[l) 

<5l,l) = (Ti,i, <5i,2> = ^(CTl.l +Cri,2)tanh2fc, <5i,3> = CTi,2, (2.11) 

and ^ 

(■si,isi,2> = -(l + o-i,io-i,2)tanh2fc, <5i, 1X1,3) = o-i,io-i,2 . (2.12) 

For the spins from the white rows 

<S2j) = 0, {s2,jS2,j+n) = tanh" k. (2.13) 

Now it is relatively easy to find the renormahzed Hamiltonian J^'Ccr) (2.7) in the form 

^'((j) = Giki) + k'Y^Gi^jGi^j+i + fc;^(7,-jo-,-+ij +n((T), (2.14) 

where G(fc, ) is a constant (independent of effective spins tj) term which can be used to calculate the free 
energy per site according to the formula 

- Gikf^) 

f-Y.—^' (2.15) 

and "n" numbers the LPRG steps. 0((7) denotes the additional interactions generated eventually by LPRG 
transformation. 

The approach presented above can be also used to consider a quantum spin model or an interacting 
electron model. However, in these cases the LPRG does not start with an exact but with an approximate 
decimation for one-dimensional systems We should also emphasize that generally, the approach is 
relevant for higher temperatures. Using LPRG one can show the existence of the critical point, and can 
obtain the location of a transition point, if any, and the temperature or field dependences of the thermo- 
dynamic quantities. 



3. Extended ANNNI model in a field 

There exists a class of compounds with a layered structure and strong c-axial anisotropy which can 
be described by the Ising-like ANNNI model or its extensions. For example, models of this kind have 
been proposed to describe the phase diagrams of the uranium compounds: UAsi-xSe^; with x < 0.1 fl3l 
or UT2Si2 (T - Pd, Ni) (l^ES made of ferromagnetic layers which comprise three ordered phases: anti- 
ferromagnetic (H 1 — ), ferrimagnetic (+ + -) and incommensurate. Each of this system exhibits a triple 
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point in which two ordered phases (ferrimagnetic and incommensurate) meet the paramagnetic one, and 
all undergo the high field phase transition from paramagnetic to ferrimagnetic (+ + -) phase. However, 
if in UAsi_j;Sex the para-ferri phase transition is of the first order, in UT2Si2 this transition, at least in 
the vicinity of the multicritical point, seems to be a continuous one. For the symmetry reason, the para- 
ferrimagnetic (+ + -) phase transition in the presence of an external field should be discontinuous and a 
question has arised if the ANNNI model can exhibit in such a case a critical transition. Within the MFA it 
has been shown that the ANNNI model should be supplemented by some additional, e.g., four-spin inter- 
action to exhibit an isolated critical point [itIi . In order to answer this question, we have used the LPRG 
to study the extended ANNNI model in a field 



('7> di) i'i) 

+ fc4 X Si,jSi+ijSij+iSi+ij+i + h^Si. 



(3.1) 



Hi) 



In the ground state in zero field, the standard ANNNI model with A;4 = 0, fci , ^2 < and jc = A;2 / fci = 0.5 

exhibits a multicritical point where antiferromagnetic phases (H h-) and (+ H ) meet a ferrimagnetic 

phase with the spins of the ferromagnetic planes ordered in the sequence (+ + -) (<2, 1) phase in the 
notation of the paper (3). The external magnetic field extends the (+ + -) phase in the {k,H) plane 
flill . All magnetic structures of UPd2Si2 as well as UNi2Si2 and UAsi_j;Sej; (0 ^ x ^ 0.1) in the fields along 
c axis exhibit ferromagnetic layers with moments perpendicular to the layers. So, only the magnetic 
order between the ferromagnetic layers is of interest. Consequently, it seems that to understand the main 
feature of the paramagnetic-ferrimagnetic (+ + -) phase transition one can confine oneself to consider the 
two dimensional (2D) problem. Thus, in the present paper we use the LPRG method to study the extended 
ANNNI model in two dimensions. 

The idea of the LPRG for ANNNI-type model in 2D with the following projector for one decimated row 



P(0-, S) = i (1 + Si,3) (1 -H CTi,2Si,6) 
4 



(3.2) 



is presented in figure [2] The full circles represent the spins which survive in the decimation procedure. 
According to figure [2] in each step of the RG transformation every other row ("even row") is removed, 
and from odd rows every third spin survives. The cluster presented in figure[2]can be used to study sys- 
tems made of ferromagnetic or antiferromagnetic chains but it does not preserve the antiferromagnetic 
and much less modulated order between the chains. However, this cluster preserves the ferrimagnetic 
(+ + -) ordering which is a matter of interest in this paper. In order to take into account both possibilities 
(antiferromagnetic and ferrimagnetic orders) one should consider a cluster of 11 chains and for three 
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Figure 2. Cluster (8-10-8-10-8) used to get renormalized HamUtonian of the extended ANNNI model in the 
LPRG procedure. 
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possibilities (antiferromagnetic + - +-, antiferromagnetic + + — , and ferrimagnetic orders) 16 chains. 
Although calculation with such clusters is straightforward it becomes quite involved. It should be em- 
phasized once more that the LPRG with the cluster presented in figure[2]is suitable to describe, except for 
para-ferro, only para-ferrimagnetic (+ + -) phase transitions. Such phase transitions can occur in both 2D 
and 3D systems. Thus, we have confined ourselves to consider 2D problem although the incommensurate 
phases observed in the mentioned uranium compounds can be described by 3D ANNNI model. 

As usual, the RG procedure leads from the set of original parameters of the Hamiltonian d2.3t (h, k, 
fci, kz, ^4) to the set of renormalized parameters (/«', k', k[, k!^, k'^) 

h'u), k'afuf, k[u{a^i^'\ k'^afa^i^^\ k'.af a^l^'^ af G'-i^'\ (3.3) 

However, in the lowest nontrivial order of the cumulant expansion, which in our case is the second order, 
several new odd and even interactions come into play. For the sake of simplicity, we will take into account 
the contributions up to the second order only from the one- and two-spin interactions, contributions to 
the first order from three- and four-spin interactions, and neglect any contributions from the higher order 
spin terms (ill]. Accordingly, we consider 17 parameters, i.e., five original ( 13. 3t and 12 generated by the 
RG transformation. 

To evaluate the transformation ( 12. 8> one has to know the averages of spins (S,-,y> and two spins prod- 
ucts {Sj^j Sij+n) from the decimated ("odd") and removed ("even") rows. The chain averages for the spins 
from decimated rows of the Ising model in a field with the weight operator <3.2) have been presented in 
pap er Q], and for the spins from the removed rows, these averages are known exactly (see for example 
11811 and references therein). Now, we are able to numerically evaluate the renormalization transforma- 
tion ( 12. 8> which in our approximation has a form of 17 recursion relations. As usual, in order to deter- 
mine the critical temperature, one has to find a critical surface which separates in the parameter space 
the region of attraction of the two stable fixed points, zero temperature, fc^ = 00 and infinite tempera- 
ture ka = 0. The existence of such a surface means that the system can undergo a second order phase 
transition. Moreover, we can also numerically calculate the free energy per spin collecting the constant 
terms generated in each step of the iteration process ( I2.15> . In a disordered phase we obtain a rapidly 
convergent infinite series for the free energy which can be used to find the specific heat as a function of 
the reduced temperature f = T/J. 

Let us start with the standard ANNNI model with fc = 1, fci = -0.6, fca = -0.3 and A;4 = The value 
of the critical temperature has been determined on the ground of the recursion relation analysis. Using 
the formula <2.15K the free energy per site and the specific heat as functions of temperature and a field 
have been found. In figure [S] the temperature dependence of the specific heat for several values of the 
field is presented. The specific heat divergence corresponding to a critical point is visible for h = and for 
h around 0.11. For < h < he ~ O.ll, the specific heat exhibits a maximum. Unfortunately, the LPRG fails 
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Figure 3. (Color online) The specific lieat temperature dependences of tlie ANNNI model witli ki = -0.6, 
^2 = -0.3, ^4 = for several values of the field. 
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Figure 4. (Color online) The specific heat temperature dependences of the extended ANNNI model with 
ki = -0.6, k2 = -0.3 at zero field for several values of the four spin coupling k^. 



2000 r 



1500 



500 - 



I — 
1.55 




1.60 



1.65 



Figure 5. (Color online) The specific heat temperature dependences of the extended ANNNI model with 
ki = -0.6, k2 = -0.3, ^4 =0.2 for several values of the field. 
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Figure 6. (Color online) Comparison of the specific heat curves of the standard ANNNI model at fields 
ft = and ft = 0.11 around the critical field with such a curve for the model with k^ = 0.2 around the 
critical field ft = 0.121. 
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Figure 7. (Color online) The temperature dependences of the nearest-neighbor <x, s,+i> (bottom curve) 
and next-nearest-neighbor ■5,-1-2) chains correlation functions for extended ANNNI (ki = 0.2) model at 
/z = 0.121. 

to properly describe the specific heat behavior below the maximum because it is not capable of describing 
the ordered phase which is expected below the maximum, if the system undergoes a discontinuous phase 
transition. For some value of the field h > he (h - 0.15 in figure[3), specific heat has only a broad hump. 

This means that the considered ANNNI in the presence of the external field can exhibit a continuous 
phase transition to the ferrimagnetic phase (-1- -1- -) at ^ = and h - he ~ 0.11. Thus, within the LPRG 
unlike in MFA it is not necessary to supplement the original ANNNI model by an additional four-spin 
interaction to reach an isolated critical point for a finite field around h-0.11 12 ill . 

Now, we proceed to the extended ANNNI model with A; = 1, fcl = -0.6, k2 - -0.3 and + 0. In figure|4] 
we present the temperature dependences of the specific heat in zero field for several values of negative 
and positive values of ^4. As seen, the four-spin interaction simply shifts the critical point towards higher 
temperature for ^4 > and towards lower temperature for ^4 < 0. Figure [s] shows the temperature de- 
pendence of the specific heat for the model with ka, - 0.2 in finite fields. Similarly to the standard ANNNI 
model case, the application of the external field changes the divergence of the specific heat into a maxi- 
mum for < h < he. For h = 0.1214, the divergence appears again indicating a critical point. In figure[6] 
we compare the shapes of the specific heat curves for the field strengths close to the critical values for 
standard and extended ANNNI models. In the case of the k4 - 0.2 model, a small peak in specific heat 
appears above the critical temperature. In figure [7] the two-site correlation functions for the spin from 
adjacent (Gi = (SiSj+i)) and next nearest neighbor (G2 - (SiS,+2>) chains are presented. As seen, a kink 
of G2 around the temperature of the specific maximum is observed. It can suggest some preliminary or- 
dering between the spins from the next nearest neighbour chains. To summarize, it has been shown that 
within the LPRG approximation, both the original ANNNI model and the one supplemented by the four- 
spin interaction with k2/ki - 0.5 and fci, fc2 < undergo the continuous phase transition at h-0 and can 
undergo such a transition for some finite value of the field h - he- In the calculations presented in this 
paper, a sharp anomaly in specific heat is evident at fields around he which could explain the experimen- 
tally observed anomalies in UPd2Si2. The four-spin interaction shifts both the critical temperature and 
field, and changes the shape of the specific heat curves. 
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MeTOA-niHiMHoY nepTyp6aTi/iBHoY peHopMa/iisaL^mHoY rpyni/i 
f\na i3i/iHronoAi6Hi/ix cniHOBi/ix ci/ictcm 

l7l. LLlHaMfl 

iHCTWTyi HH3bK0i TeMnepaiypn i crpyKTypHnx flOC/iiflxeHb, Ho/ibCbKa aKafleMm HayK, Bpoi^/iaB, no;ibLL(a 

/lihiPiHe nepTyp6aTWBHe nepeiBopeHHS (LPRG) BUKopwcTOByeibca fl/ia BWB4eHHfl lepMOflWHaMiKn aKcia/ibHOi 
MOfleyii iBWHra 3 HacTynHi/iMW flo Haw6;inx4nx cyciflaMU 3 HOTupwcniHOBOio B3aeM0flieK) {po3LunpeHa MOfle/ib 
ANNNI) y noyii. HpeflCTaB/ieHO LPRG fl/ia c;ia6oB3aeMOfliK)4nx ;iaHL4K))KKiB l3i/iHra. Merofl 3acT0C0BaH0 flo bh- 
B4eHHa napa-c|)epoMarHiTHHx c(>a30Bnx nepexofliB y CKiH^eHHOMy no.ni, mo cnocrepiraiOTbCH b LuapyBarux cno- 
yiyKax ypany UAsi_;cSe;c, UPd2Si2 4H UNi2Si2. Bume sraflani cwcreMn 3po6;ieHi 3 c|)epoMarHiTHHx luapiB i cni- 
Hi/i 3 HaPi6;inx4nx i HaciynHwx flo Hati6;in>K4MX LuapiB e 3B'fl3aHi aHTUcfiepoMarHiTHUMM B3aeM0fli;iMH /i < 
i /2 < 0, BiflnoBiflHO. Ko>KHa 3 L^nx cucreM fleivioHCTpye norpiPiHy T04Ky, b m\\A p,B\ Bnop^flKOBaHi c|)a3n (c|)e- 
poMarHiiHa i HecniBMipHa) 3ycTpiHaK)TbC5i 3 napaMarHiTHOK) c(>a30Kj i Bci c|)a3n 3a3HaiOTb (|)a30Boro ^epexo/^y 
y ci/i;ibH0My no/ii 3 napa- flo c|)epoMarHiTHOi' (-1- -I- -) (t)a3w. flpoTe, flKmo b UAsi-xSex £ napa-c|)epo (|3a30BwCi 
nepexifl nepiuoro pofly, OHiKyeibcn 3 cuMerpiRHHX MipKyBaHb, b UT2Si2 (T = Pd, Ni) u,ei/\ nepexifl BHflaeTbC5i 
HenepepBHUM, npuHawMHi b OKO/ii My;ibTMKpnTi/i4H0TT04KH. B paMKax Ha6;inxeHHfl cepeflHboro nonsf, K.pmm- 
Hi/iti xapaKiep napa-c|)epoMarHiTHOro nepexofly b CKiH^eHHOMy no/ii, npuHaRMHi b OflHiPi i3o;ibOBaHiM T04i^i, 
MO>Ke 6yTH oni/icaHWM 3a flonoMoroK) MOfle/ii ANNNI, «Ka flonoBHeha flOflaiKOBOK), HanpuK/iafl, HOTupwcniHO- 
BOK) B3aeM0flieK). npore, b Ha6;inxeHHi LPRG fl/is Koe(()iL4ieHTa k = /2//1 no6/in3y 0.5 e KprnnHHe 3Ha4eHHfl 
no/ia, flyifl flKoro i3o;ibOBaHa KpnTH4Ha T04Ka laKOx iCHye b opuriHa/ibHm MOfle/ii ANNNI. no3HTHBHa 40Twpn- 
cniHOBa B3aeM0fliH 3cyBae Kpi/irnHHy roHKy flo bi/iihi/ix no/iis i 3MiHK)e (jjopwy KpuBoT nmoMoT len/ioeMHOCTi. B 
ocraHHbOMy BunaflKy fl/ia flOcrarHbO Ma/ii/ix no/iiB nmoMa len/ioeMHiCTb fleMOHCrpye flBoniKOBy crpyKiypy b 
napaMarhiTHiti <^a3\. 

KniOMOBi cyiOBa: MOflejib ANNNI, peHopMayiisaLfiUHa rpyna, iBOJibOBana Kpmi/mna roHKa 
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